What environmental signals prompt native fish to leave the floodplain before disconnection with the river and avoid stranding?
The goal of my current research is to help understand and conserve native salmon (Oncorhynchus tshawytscha) populations by better understanding fish utilization of natural, restored, or managed floodplains under varying environmental conditions. The lower Cosumnes River is an ideal place to research these dynamics, as it has no major dams, so the flow is relatively natural, and there is a restored floodplain that can be used to understand timing and use by juvenile salmon and other native fishes. Information from my research will inform managers and help prioritize restoration efforts throughout the Central Valley.
library(tidyverse)
library(plotly)
glimpse(storms)
## Observations: 10,010
## Variables: 13
## $ name <chr> "Amy", "Amy", "Amy", "Amy", "Amy", "Amy", "Amy", "...
## $ year <dbl> 1975, 1975, 1975, 1975, 1975, 1975, 1975, 1975, 19...
## $ month <dbl> 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7,...
## $ day <int> 27, 27, 27, 27, 28, 28, 28, 28, 29, 29, 29, 29, 30...
## $ hour <dbl> 0, 6, 12, 18, 0, 6, 12, 18, 0, 6, 12, 18, 0, 6, 12...
## $ lat <dbl> 27.5, 28.5, 29.5, 30.5, 31.5, 32.4, 33.3, 34.0, 34...
## $ long <dbl> -79.0, -79.0, -79.0, -79.0, -78.8, -78.7, -78.0, -...
## $ status <chr> "tropical depression", "tropical depression", "tro...
## $ category <ord> -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, ...
## $ wind <int> 25, 25, 25, 25, 25, 25, 25, 30, 35, 40, 45, 50, 50...
## $ pressure <int> 1013, 1013, 1013, 1013, 1012, 1012, 1011, 1006, 10...
## $ ts_diameter <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA...
## $ hu_diameter <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA...
library(viridis)
ggplot(storms) +
geom_jitter(aes(wind, pressure, color = status)) +
scale_color_viridis_d() +
labs(title = "Hurricane Pressure v Wind", subtitle = "By Status") +
theme_bw()
Create a table that identifies the mean wind, pressure, ts_diameter, hu_diameter of each status of storm.
library(htmlTable)
storms_tbl <- storms %>%
drop_na(ts_diameter, hu_diameter) %>%
group_by(status) %>%
summarize(mean_wind = mean(wind), mean_pressure = mean(pressure), mean_ts = mean(ts_diameter), mean_hu = mean(hu_diameter))
storms_tbl$mean_wind <- round(storms_tbl$mean_wind, 2)
storms_tbl$mean_pressure <- round(storms_tbl$mean_pressure, 2)
storms_tbl$mean_ts <- round(storms_tbl$mean_ts, 2)
storms_tbl$mean_hu <- round(storms_tbl$mean_hu, 2)
htmlTable(storms_tbl)
| status | mean_wind | mean_pressure | mean_ts | mean_hu | |
|---|---|---|---|---|---|
| 1 | hurricane | 87.15 | 966.35 | 288.11 | 72.96 |
| 2 | tropical depression | 28.21 | 1006.47 | 0 | 0 |
| 3 | tropical storm | 45.75 | 999.03 | 159.61 | 0.04 |